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ABSTRACT 

We present results from Chandra, XMM-Newton, and ROSAT observations of the Planck SZ- 
detected cluster A37I6 (PLCKG345.40-39.34 - G345). We show that G345 is, in fact, two subclusters 
separated on the sky by 400 kpc. We measure the subclusters’ gas temperatures 2-3 keV), total 
1-2 xlO®'^ Mq) and gas (~ 1”2 xlO®® Mq) masses, gas mass fraction within rsoo, entropy profiles, 
and X-ray luminosities 10^® erg s“®). Using the gas density and temperature profiles for both 
subclusters, we show that there is good (O.Scr) agreement between the expected Sunyaev-Zel’dovich 
signal predicted from the X-ray data and that measured from the Planck mission, and better agree¬ 
ment within 0.6 ct when we re-computed the Planck value assuming a two component cluster model, 
with relative amplitudes fixed based on the X-ray data. Dynamical analysis shows that the two galaxy 
subclusters are very likely (> 97% probability) gravitationally bound, and in the most likely scenario, 
the subclusters will undergo core passage in 500 ± 200 Myr. The northern subcluster is centrally 
peaked and has a low entropy core, while the southern subcluster has a high central entropy. The 
high central entropy in the southern subcluster can be explained either by the mergers of several 
groups, as suggested by the presence of five giant ellipticals or by AGN energy injection, as suggested 
by the presence of a strong radio source in one of its massive elliptical galaxies, or by a combination 
of both processes. 

Subject headings: galaxy clusters: general — cosmology: large-structure formation 


1. INTRODUGTION 

Glusters of galaxies are the largest gravitationally 
bound structures in the Universe to have reached virial 
equilibrium. In the standard ACDM cosmology, mas¬ 
sive halos dominated by dark matter assemble by ac¬ 
cretion of smaller groups and clusters. Under the 
influence of gravity, uncollapsed matter and smaller 
collapsed halos fall into larger halos and, occasion¬ 
ally, halos of comparable mass merge with one an¬ 
other. Observations of substructures in clusters of galax¬ 
ies (e.g., Jones & Forman 1984, 1999; Mohr et al. 1995; 
Buote & Tsai 1996; Jeltema et al. 2005; Lagana et al. 
2010; Andrade-Santos et al. 2012, 2013) and the growth 
of structure (Vikhlinin et al. 2009) strongly suggest 
that clusters formed recently (Forman & Jones 1982; 
Richstone et al. 1992). To better understand the evo¬ 
lution of the Universe, it is important to identify and 
characterize these large structures. 

The first catalog of 189 Sunyaev-Zel’dovich (SZ) clus¬ 
ters detected by the Planck mission was released in 
early 2011 (Planck Collaboration et al. 2011). Through 
a Chandra XVP (X-ray Visionary Program - PL Chris¬ 
tine Jones) and HRC Guaranteed Time Observations 
(PI: Stephen S. Murray), which comprise the Chandra- 
Planck Legacy Program for Massive Glusters of Galax¬ 


ies®, we are obtaining Chandra exposures sufficient to 
collect at least 10,000 source counts for each of the ESZ 
Planck clusters at z < 0.35. PLCKG345.40-39.34 (here¬ 
after G345 - RA = 20:52:16.8, DEC = -52:49:30.7) is a 
nearby cluster (z = 0.045). G345 was the first cluster ob¬ 
served as part of the Chandra XVP, which revealed it to 
be a double cluster (we distinguish the north and south 
subclusters as G345N and G345S, respectively). 

Historically, Abell et al. (1989) had catalogued a sin¬ 
gle cluster of galaxies (A3716 - cross mark in Figure 1) 
located ~ 6' northwest of G345S and ^ 6' southwest of 
G345N. 

In Section 2 of this paper we present the Chandra, 
XMM-Newton, and ROSAT observations and data re¬ 
duction. In Section 3 we describe the X-ray spatial and 
spectral analysis. The total mass and gas mass for each 
subcluster are presented in Section 4, and the entropy 
profiles of both subclusters are computed in Section 5. In 
Section 6 we compute the expected Sunyaev-Zel’dovich 
signal from the X-ray data, as well as a comparison be¬ 
tween the X-ray emission and the Planck reconstructed 
Sunyaev-Zel’dovich map. Finally, we present a dynam¬ 
ical model for the two subclusters in Section 7 and the 
conclusions in Section 8. 

1 hea-www.cfa.harvard.edu/CHANDRAJLANCK.CLUSTERS/ 
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Fig. 1.— Left: 0.5-2.0 keV XMM-fVeuiton image of G345. Right: ESO DSS-Red optical image of the same sky region 
presented in the left panel. Isointensity X-ray contours are shown in both panels. The cross indicates the Abell et al. (1989) 
catalogue position for A3716. 


We assume a standard ACDM cosmology with LIm = 
0.3, IIa = 0.7 and i7o = 70 km s“^Mpc“ , implying a 
linear scale of 0.89 kpc arcsec“^ at the G345 luminosity 
distance of 199 Mpc (z = 0.045). 

2. X-RAY OBSERVATIONS AND DATA 
REDUCTION 

2.1. Chandra 

We observed C345 on November 21 and 26, 2012, with 
the Chandra X-ray Observatory (ACIS-I detectors, VF 
mode, Obsids 15133 - 15 ks - and 15583 - 15 ks - see 
center panel of Figure 2). The data were reduced follow¬ 
ing the processing described in Vikhlinin et al. (2005), 
applying the calibration files CALDB 4.5.3. The data 
processing includes corrections for the time dependence 
of the charge transfer inefficiency and gain, and a check 
for periods of high background (none were found). Also, 
readout artifacts were subtracted and standard blank sky 
background files were used for background subtraction. 

2.2. AMM-Newton 

We observed C345 on April 2, 2013, with the XMM- 
Newton Observatory for 22.7 ks (Obsid 0692930101). 
We used the images script^ from the XMM-Newton web¬ 
site to create exposure corrected images of the cluster in 
the 0.5-2.0 keV energy band. This script removes periods 
of high background, as well as bad pixels and columns, 
spatially smooths and exposure corrects the data and 
merges PN and MOS observations (see Figure 1 left panel 
and Figure 2 right panel). 


^ xmm.esac.esa.int/external/xmm_science/gallery/utils/images.shtml 


2.3. ROSAT 

C345 was observed on April 2, 1991 for 2.7 ks (Obsid 
RP800042A00) and on April 1, 1992 for 1.6 ks (Obsid 
RP800042A01) with the ROSAT PSPC. We used avail¬ 
able hard band (0.4-2.4 keV) images, background, and 
exposure maps to create a background subtracted, ex¬ 
posure corrected image of the cluster that we used to 
constrain the subclusters’ surface brightness profiles at 
large radii (see Figure 2 left panel). 


3. SPATIAL AND SPECTRAL ANALYSIS 
3.1. X-ray Surface Brightness Radial Profiles 

The surface brightness is the projection of the plasma 
emissivity along the line of sight. We fit the X-ray surface 
brightness radial profile of each C345 subcluster with a 
/3-model (Cavaliere & Fusco-Femiano 1976) which is well 
suited for non-cool-core relaxed clusters. This model is 
defined as: 


S{r) = So 


1 + ( - 

rn 


-3/3-1-0.5 


( 1 ) 


where r^ is the core radius, /3 is the shape parame¬ 
ter, and So is the central surface brightness. Figure 3 
shows the surface brightness profiles for the northern 
(left panel) and southern (right panel) subclusters, along 
with the best /3-model fits. We detect X-ray emission 
from Chandra and XMM-Newton to a radius of ~ 700 
kpc for the northern component and to ~ 1 Mpc for the 
southern one. The background level is (1.70 ± 0.13) x 
10“^ counts s“^ arcsec“^ (normalized to the Chandra 
ACIS-I count rate sensitivity in the 0.5-2.0 keV energy 
band). For C345N, we modeled the X-ray surface bright¬ 
ness with the addition of a Caussian function to describe 
the X-ray emission associated with the central giant el- 
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TABLE 1 

Best Fit Parameters for the /3-Model and central density. 


Subcluster 

rc 

0 

n-efi 

PO 


(kpc) 


(cm“®) 

(g cm-3) 

G345N 

118 ± 9 

0.54 ± 0.02 

(3.19 ±0.28) X 10-3 

(6.21 ±0.54) X 10-2’^ 

G345S 

228 ± 12 

0.48 ± 0.01 

(8.96 ±0.51) X lO-'^ 

(1.74 ±0.10) X 10-2’^ 


Note. — Columns list subclusters’ names, determination and Icr uncertainties for the core radius, /3 (Equation (1)), and central electron 
number and gas mass densities (Equation (4)). 


TABLE 2 

Comparison Between the Best Fit Parameters for the /3-Model Using Different Data Sets. 


Telescope 

Subcluster 

(kpc) 

0 

X?ed 

Chandra 

G345N 

82 ± 23 

0.44 ±0.08 

0.83 


G345S 

228 ± 49 

0.44 ± 0.04 

1.37 

XMM- Newton 

G345N 

91 ±33 

0.46 ±0.10 

1.05 


G345S 

199 ± 20 

0.45 ±0.02 

1.73 

Chandra + X^MNL-Newton + ROSAT 

G345N 

118 ± 9 

0.54 ± 0.02 

1.07 


G345S 

228 ± 12 

0.48 ± 0.01 

1.31 


Note. — Columns list telescope data used for fitting the surface brightness profile, subclusters names, best fit for the core radius, /3 
(Equation (1)), and reduced x-square. 



Fig. 2.— Left: The 0.4-2.4 keV ROSAT PSPC image of G345. Center: The 0.5-2.0 keV Chandra image of G345. Right: The 
0.5-2.0 keV XMM-Newton image of G345 showing rsoo for each subcluster. All three images have the same spatial scale. 


liptical (ESO 187-G026). The Gaussian function in the 
Chandra Ht describes the bright core, whereas in the case 
of jCMM-N ewton fit, the Gaussian function is broader 
due to the telescope’s larger PSF compared to Chandra. 
With Chandra, we extracted a spectrum from a circle 
of 7 kpc (7.9") in radius centered on ESO 187-G026 to 
estimate the gas temperature of the bright central re¬ 
gion (using the spectrum from an annulus with inner 
and outer radii of 7 (7.9") and 27 kpc (30.5") respec¬ 
tively, as the background component to properly sub¬ 
tract the cluster emission). Eor the galaxy, we obtain 
kT = 0.92toll keV and Lx = x erg s-b In 

the Appendix we show that the (stellar -|- LMXB) lumi¬ 
nosity of this galaxy is Lq. 5 - 2.0 keV = l-fd x 10^^° erg s“^ 
which is only ~ 7% of the total 0.5-2.0 keV luminos¬ 
ity within 7 kpc (7.9") of its center. Thus, the major¬ 
ity of the X-ray emission comes from diffuse gas. In a 
similar way, Vikhlinin et al. (2001) showed that there 


are ^ 3 kpc cores of 1-2 keV gas in both giant ellip¬ 
ticals (NGC4874 and NGC4889) in the massive and hot 
{kT ^ 8 keV) Goma cluster, and Sun et al. (2007) pre¬ 
sented a systematic investigation of X-ray thermal coro- 
nae in a survey of 25 nearby clusters. They showed that 
cool galactic coronae {kT = 0.5-1.0 keV generally) are 
common in cluster cores. Using ROSAT’s larger field of 
view to better constrain the background level along with 
Chandra and XMM-Newton observations, we fitted each 
subcluster X-ray surface brightness to a /3-model pro¬ 
file and determined /3 and Tc (the background level for 
each surface brightness profile {Chandra, XMM-Newton, 
and ROSAT) was free to vary independently, while r^ 
and /3 were tied to a single best fit value). F^or G345N, 
we obtained /3 = 0.54 ± 0.02 and Tc = 118 ± 9 kpc 
(Xred = ^'^d for G345S, we obtained /3 = 0.48 ± 

0.01 and Tc = 228 ± 12 kpc (Xred = 1.31). We also 
fit the Chandra and XMM-Newton profiles individually 
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and found consistent /3 and re- The best fit parameters 
are summarized in Table 1. The southern subcluster has 
a flatter (3 and a larger core radius than the northern 
subcluster, which along with the poorer /3-model fit and 
higher central entropy of G345S (the entropy profiles for 
each subcluster will be discussed in Section 5) suggest 
this subcluster is still in the process of forming. How¬ 
ever, visual inspection of Figure 3 as well as the best 
fit reduced show that the /3-model describes well the 
X-ray surface brightness profiles of G345N (with the ex¬ 
ception of the central region due to the central galaxy’s 
cool galactic corona X-ray emission) and G345S, validat¬ 
ing the assumption that the gas density defined by the 
surface brightness follows this model. 


3.2. Gas Density Radial Profiles 

The /3-model gas density distribution that corresponds 
to the surface brightness distribution given by Equation 
( 1 ) is: 

-3/3/2 

( 2 ) 


(f) = Pg,0 


l-f 


where Pg^o is the central gas density. The core radius, 
Tc, and /5 are constrained from fitting the X-ray surface 
brightness profile. 

The gas mass within rsoo (the radius defining a sphere 
whose interior mean mass density is 500 times the critical 
density at the cluster redshift - see Section 4) is then 
given by: 


Afg,500 — 47rpg_o 


47r 


prsoo 

f r\ 

/ 

1 + (- I 

Jo 

\rcj 


21 -3/3/2 


^dr 


= fPa,orloo 2i"i(3/2, 3/3/2; 5/2; -y^), (3) 


where y = and 2 ^ 1 ( 0 , 6 ; c; d) is the Gauss hyper¬ 

geometric function. For a plasma with a given electron 
to hydrogen number densities ratio rie/uH, the central 
electron number density, rie.o is calculated as: 


'^e,0 — 


rK{nJnH)[DA{l + z)]H6P - 3) x lO^^x 


1 -b 


26 / 3 - 3 ^ 35 ( 3 ^ - 1/2,3/3 - 1/2) 

-3/3+3/2 |- 2l -3/3-I-3/2\ 

l+(^l I (4) 


where Da is the angular distance of the cluster, B{a, b) 
a, 5 G K is the Euler /3 function, i?i and Rf are the pro¬ 
jected initial and final radii of the annulus for which the 
spectrum has been ht and the normalization, K, was 
computed using the APEC model in XSPEG. The relation 
between the gas density and electron number density is 
given by pg = /ie^ieWa, where pe is the mean molecular 
weight per electron and ma is the atomic mass unit. For 
a metallicity of 0.3 Zq, using the reference values from 
Anders & Grevesse (1989) we obtain pe = 1.1706 and 
nejnn = 1.1995. The central electron number and gas 
densities for G345N are approximately four times higher 
than those for G345S. These are given in Table 1. 


3.3. Gas Temperature Radial Profiles 


The temperature profiles of most clusters have a broad 
peak within 0 . 1 - 0.2 r 2 oo and decreases at larger radii, 
reaching approximately 50% of the peak value near 
0 . 5 r 2 oo (Vikhlinin et al. 2006). In cool-core clusters 
there is also a temperature decline toward the cluster 
center due to radiative cooling. The analytic model con¬ 
structed by Vikhlinin et al. (2006) for the 3D tempera¬ 
ture profile describes these general features. However, 
due to the quality of our data, we employed a simpli- 
hed form of this temperature profile with some of the 
parameters fixed (the universal temperature profile from 
Vikhlinin et al. (2006), but with the transition radius, rt, 
allowed to vary). Thus 


T3D(r) 


To _ [r/(0.075rt)]^® + 0.45 
[1 + (r/rt)2] ^ [r/(0.075rt)]i-9-hl 


(5) 


Monte-Carlo simulations were performed to estimate 
the uncertainties in the best fit values for the parameters 
in this analytical model. This analytic model for T 3 D('r) 
(Equation (5)), allows very steep temperature gradients. 
In some realizations, such profiles are formally consistent 
with the observed projected temperatures because pro¬ 
jection flattens steep gradients. However, steep values of 
dT/dr often lead to unphysical mass estimates, for exam¬ 
ple, profiles with negative dark matter density at some 
radii. We addressed this problem in the Monte-Carlo 
simulations by accepting only those realizations in which 
the best-ht T{r) leads to ptot > Pgas in the radial range 
r < l.brsoo. Finally, in the same radial range, we ver¬ 
ified that the temperature profiles corresponding to the 
mass uncertainty interval are all convectively stable^, i.e., 
dlnT/dlnpg < 2/3. 

Using XMM-Aewton observations, we extracted spec¬ 
tra in three half annuli subtending 180 degrees on the side 
opposite to the companion subcluster, centered on each 
subcluster in the radial range from from 0 to ^ 400-600 
kpc for each subcluster. Each half annulus has at least 
3000 source counts. We fit these with an absorbed APEC 
model. The measured Galactic hydrogen column density 
in the direction of the cluster is /Vh = 2.51 x 10^° cm“^, 
which was kept fixed while fitting the spectra. We then 
followed the procedures described below to obtain the 
2D and 3D temperature profiles. The measured 2D, fit¬ 
ted 2D, and 3D temperature profiles are presented in the 
left panel of Figure 4. The 2D temperature profile was 
computed by projecting the 3D temperature weighted by 
gas density using the spectroscopic-like temperature (a 
formula for the temperature which matches the spectro¬ 
scopically measured temperature within a few percent, 
Mazzotta et al. 2004): 


^ f PlT^dz 

f Pl'^3~D'''^dz 


( 6 ) 


The spectroscopic-like temperature was also computed 
in the (0.15 - 1) rsoo range using the model for the 3D 


3 Assuming the motion of a gas element in a medium in hydro- 
static equilibrium is adiabatic, one can apply Newton’s second law 
to the net force per unit volume of the gas to obtain two solutions: 
an oscillatory one, which is stable, and a run-away motion, which is 
unstable. The stable solution imposes dlnT/dlnpg <7 — !, where 
7 is the adiabatic index. This is also known as the Schwarzchild 
criterion. For 7 = 5/3, dlnT/dlnpg < 2/3. 
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Fig. 3.— Left: X-ray background-subtracted surface brightness profiles and fits for G345N (solid lines), along with the surface 
brightness fit for G345S (dashed line). The difference between the profiles in the core is due to the XMM-Newton poorer angular 
resolution compared to that of Chandra. The Gaussian function in the Chandra fit describes the cool galactic corona. For 
XMM-iVeuiton the Gaussian function to describe the core is broader, corresponding to the telescope’s larger PSF. Right: X-ray 
surface brightness profile of G345S (dashed line), along with the surface brightness fit for the G345N cluster component (solid 
line) for comparison. For each profile, regions in the direction of the other subcluster are excluded. For display purpose, the 
XMM-A^euiton central surface brightness in the S-model was normalized to match the central surface brightness in the Chandra 
fit. The background level is (1.70 ± 0.13) x 10“~ counts s“^ arcsec”^ (normalized to the Chandra ACIS-I count rate sensitivity 
in the 0.5-2.0 keV energy band). 


temperature and the gas density profile. For G345N, we 
found fcr= 2.11±0.16keV, andforG345S, fcT= 2.85± 
0.20 keV. 


4. TOTAL AND GAS MASSES OF SUBGLUSTERS 

Given the three-dimensional models for the gas density 
and temperature profiles, the total mass within a radius 
r can be computed for each subcluster from the equation 
for hydrostatic equilibrium (e.g., Sarazin 1988), 


M(r) = 


—kTr / dlnpg 


^itihG \ d\nr 


= -3.67 X 10^^ MQkTr 


dlnT\ 
d\nr J 
f din pg 


\ dinr 


dlnT\ 
dlnr / 


(7) 


where k is the Boltzmann constant, T is the gas tem¬ 
perature in units of K, and r is in units of Mpc (the 
normalization corresponds to p = 0.6107, appropriate 
for 0.3 Zq). 

Using Equation (7), r^oo is computed by solving 


we compute a gas mass within rsoo of Mg^soo = (1.45 ± 
0.23) X IO^^Mq for G345N, and Mg ,500 = (1.45 ±0.14) x 
IO^^Mq for G345S. 

The gas and total masses and density prohles, as well 
as the gas mass fractions within r^oQ for both G345N and 
G345S are shown in the center and right panels of Figure 
4. The gas mass fraction within rsoo is /g = 0.084±0.008 
for G345S and is in agreement with the expected value 
from Vikhlinin et al. (2006) for clusters ~ 1-2 x lO^^M© 
(fg ^ 0.09), whereas the gas mass fraction for G345N 
is slightly larger, fg =0.130 ± 0.023. We measure the 
X-ray luminosity within 245 kpc and extrapolate to t^qq 
for both subclusters, using the best fit /3-model parame¬ 
ters. For G345N, we obtain a bolometric X-ray luminos¬ 
ity Lx.r^oo = (2.65 ±0.37) x 10^^ erg s“^, and for G345S, 
we obtain Lx,moo — (3.72 ± 0.51) x 10^^ erg s“^. These 
results are summarized in Table 5. 

Alternatively, if we use the gas mass and temperature, 
the total mass can be estimated from the Yx-M scaling 
relation of Vikhlinin et al. (2009), 


M(r5oo) = 500p,(4^/3)r|oo, (8) 

where pc is the critical density of the Universe at the 
cluster redshift. For G345N, we obtain r^oo = 724 ± 
23 kpc and a corresponding hydrostatic mass of 
Msoo.hyd = (1-13 ± 0.11) x IO^^Mq. For G345S, we 
find rsoo = 836 ± 25 kpc with a corresponding hydro¬ 
static mass of Msoo.hyd = (1.73±O.15 )x 1O14M0. Using 
the best ht parameters for the /3-model (see Table 1), 


( Y \ Bym 

3xlOi40kev) ’ 

(9) 

where Yx = Mg ,500 x fcTx, Mg ,500 is given by Equa¬ 
tion (3), and Tx is the spectroscopic-like temperature 
(see Equation ( 6 )) in the (0.15-1) x rsoo range, com¬ 
puted based on the model for the gas density and 3D 
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Fig. 4.— Top: G345N. Bottom: G345S. Left: Temperature profiles. Measured projected temperatures are shown by circles with error bars. 
Solid red and blue lines show the best fit 3D model and the corresponding projected profile, respectively. Center: Mass and density profiles 
(M(< r) increases with radius while p(r) decreases). Red and blue lines show the results for the total mass and gas mass, respectively. 
Right: Gas mass fraction as a function of radius. Solid line shows the enclosed /gas = -^g(< '^)/-^tot(< The vertical line shows the 
radius rsoo derived from the best fit mass model assuming hydrostatic equilibrium. The dashed lines show the 68% confidence level limits 
for the best fits. 


TABLE 3 

Temperature Determination 


Subcluster 

r (kpc) 

kT (keV) 

Subcluster 

r (kpc) 

kT (keV) 


14.5 - 58.5 

2.01 ±0.15 


0 - 106.8 

3.39 ± 0.25 

G345N 

58.5 - 172.7 

2.15 ± 0.11 

G345S 

106.8 - 186.0 

3.71 ± 0.31 


172.7 - 591.5 

2.21 ± 0.27 


186.0 - 398.4 

3.08 ± 0.26 


Note. — Columns give subclusters’ names, radial range for the annuli used for temperature extraction, and the measured temperatures 
with 1(7 uncertainties. 


TABLE 4 

Parameters for the Temperature Profile 


Subcluster 

To 

n 


(keV) 

(kpc) 

G345N 

2.74 ± 0.14 

1266 ± 243 

G345S 

4.07 ± 0.22 

1459 ± 141 


Note. — Columns list best fit values for Tq and rt given by 
Equation (5). 


temperature profiles. Aym = 5.77 x and 

Bym = 0.57 (Maughan et al. 2012). Here, Myx.soo is 
the total mass within rsoo, and E{z) = [Hm(1 + + 

(1 — Hm — 11a)(1 + z)'^ + is the function describ¬ 


ing the redshift evolution of the Hubble parameter. As 
for the mass determination assuming hydrostatic equi¬ 
librium, we estimated Icr uncertainties in the lx derived 
quantities using Monte Carlo simulations. We also added 
to the Monte Carlo procedure a Icr systematic uncer¬ 
tainty of 9% in the mass determination, as discussed by 
Vikhlinin et al. (2009). 

Using Equations (8) and (9), we compute rsoo.Yx = 
771 ± 47 kpc for G345N with a corresponding total 
mass of Msoo.Yx =(1.37 ± 0.25) x IO^'^Mq and gas mass 
A/g,5oo.Yx = (1.62±0.38) X IO^^Mq. For G345S, we ob¬ 
tain rsoo.Yx = 793 ± 42 kpc with a corresponding total 
mass of Msoo.Yx = (1.48±0.24) x IO^^Mq, and gas mass 
A7g,5oo,Yx = (1-33 ± 0.21) X IO^^Mq, in agreement with 
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results obtained assuming the cluster to be in hydrostatic 
equilibrium. 

Using the Yx relation, the gas mass fractions within 
?' 50 oare/g= 0.117±0.015 for G345N and/g = 0.090± 
0.007 for G345S. These values are consistent with gas 
fractions computed above assuming hydrostatic equilib¬ 
rium. These results are summarized in Table 6 . 


5. SUBCLUSTER ENTROPY PROFILES 
The entropy index of the intracluster gas is defined as 


K = 


kT 

2/3 ’ 
Ue 


( 10 ) 


where k is the Boltzmann constant, T is the gas temper¬ 
ature, and rie is the electron density. The entropy profile 
reflects the thermodynamic history of the cluster. The 
entropy increases when heat energy is deposited into the 
ICM, and decreases when radiative cooling carries heat 
energy away (Voit et al. 2005). 

To better understand the thermodynamic history of 
the G345 subclusters, we computed their entropy pro¬ 
files, which are presented in Figure 5. They are remark¬ 
ably different. G345S has extremely high entropy in its 
core (~ 200 keV cm^) due to the extremely low gas den¬ 
sity for its observed ^ 3 keV temperature (typically, the 
central entropy of a relaxed cluster with ~ 3 keV gas tem¬ 
perature would be ^ 10“30 keV cm^, Voit et al. 2005; 
McDonald et al. 2013). 

Here, we present two heating mechanisms that could 
explain an entropy increase in the ICM. In the first, clus¬ 
ters that have experienced recent major mergers can have 
high central entropies due to energy deposited in the ICM 
through shocks during the merger. The shock heating 
will produce an entropy increase, in the case of a merger 
of clusters of similar masses, that can be directly related 
to the Mach number by: 
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where Ki and K{ are the initial and final entropies, re¬ 
spectively (Zel’dovich & Raizer 1967). 

Assuming that the central entropy index of G345S was 
similar to that of G345N before it was perturbed, one 
can put a well-defined requirement on the shock strength 
in the core to boost the central entropy by a factor 
of 3.47 (ratio of central entropies between the southern 
and northern subclusters). Inserting this entropy ratio 
into Equation (11), one obtains a high Mach number of 
M = 4.8, which corresponds to ^ 4300 km s“^ for a 
3 keV cluster. This result is implausibly high (assum¬ 
ing that the mass density profile of G345S is described 
by a NEW profile with concentration parameter of 4, a 
point mass falling from infinity would arrive in the cen¬ 
ter of G345S with a velocity of ^ 2600 km s“^). How¬ 
ever, shocks are not the only way to raise the entropy 
in a merger. Dissipation of turbulence also can play an 
important role. In the case in which the gas entropy en¬ 
hancement of the southern subcluster had been caused 
by a first violent encounter with the northern subclus¬ 
ter, the cool core {kT ^ 0.9 keV - as presented in Sec¬ 


tion 3.1) in the central galaxy of G345N (ESO 187-G026) 
could be the gas that survived core passage after a colli¬ 
sion with the southern subcluster. However, such a high 
velocity encounter would have highly disturbed the X- 
ray morphologies of both G345 subclusters, which is not 
observed. 

For comparison, we searched for other clusters that 
present a similar high central gas entropy {Kq > 
150 keV cm^) and low temperature {kT < 3.2 keV) in 
the work of Cavagnolo et al. (2009), and found: AlOO, 
A193, A400, A562, A2125, and ZWCL 1215. Each clus¬ 
ter shows evidence of recent merger activity. A160 {z = 
0.0447) has two giant ellipticals in the core, that proba¬ 
bly are in the process of merging. A193 {z = 0.0485) 
has a central galaxy with a triple nucleus (1C 1695, 
Seigar et al. 2003). Since the typical time scale for mul¬ 
tiple nuclei to merge into a single one is on the order of 
a Gyr (Seigar et al. 2003), observing a triple nucleus is 
a strong indication of recent merger activity. A400 {z = 
0.024) is a well studied system which presents indications 
of merger activity (Beers et al. 1992). A562 {z = 0.11) 
has a Wide Angle Tail (WAT - radio lobes which are bent 
due to ram pressure as the host galaxy moves through 
the intracluster gas, Douglass et al. 2011). This is also 
a strong indicator of merger activity. Chandra observa¬ 
tions, together with multiwavelength data, indicate that 
the A2125 complex {z = 0.2465) is probably undergoing 
major mergers (Wang et al. 2004). Finally, inspecting 
the VLA FIRST image of ZWCL 1215 {z = 0.075), we 
notice a disturbed radio morphology associated with a 
galaxy (4C±04.41) that is probably merging from the 
southwest. However this is not strong evidence of a ma¬ 
jor merger. 

Five of the six clusters with high central gas entropy 
and low temperature show strong indications of recent or 
ongoing merger activity which is likely responsible for the 
enhancement of the central gas entropy through shocks. 
This may suggest that the high central gas entropy of 
G345S also was caused by a recent merger. However, 
if this were a merger with the northern subcluster, we 
would expect to observe an increase in the central en¬ 
tropy of the northern subcluster, which is not observed, 
therefore making this scenario unlikely. 

The core of G345S hosts hve giant ellipticals (see Fig¬ 
ure 1) (unlike G345N, which has a single dominant giant 
elliptical), suggesting that we are witnessing the forma¬ 
tion of the southern subcluster as groups of galaxies are 
in the process of merging. This scenario can explain the 
high central entropy of G345S, as the merging of the 
groups heats the gas through shocks. 

As an alternative to multiple mergers of groups, less vi¬ 
olent and frequent energy deposition into the ICM could 
also produce a significant entropy enhancement. Isobaric 
heating is a reasonable approximation for a less violent 
energy deposition where the gas pressure remains roughly 
constant while part of the energy is converted into pdV 
work, causing expansion of the gas, and part is trans¬ 
ferred to internal energy, heating the gas. For isobaric 
heating, the added heat <5(5 is related to the entropy in¬ 
crease by 
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TABLE 5 

Physical Properties Derived Assuming Hydrostatic Equilibrium 


Subcluster 

rsOO.hyd 

-^g,500,hyd 

TfsOO.hyd 

/gaSjhyd 

^Tspec 

Lx,245 

/^X.rsoo 


(kpc) 

(Mq) 

(Mq) 


(keV) 

(1043 erg s-l) 

(1043 erg s“i) 

G345N 

724 ± 23 

(1.45 ±0.23) X 1013 

(1.13 ±0.11) X 1014 

0.130 ±0.023 

2.11 ±0.16 

1.35 ±0.08 

2.65 ± 0.37 

G345S 

836 ± 25 

(1.45 ±0.14) X 10i3 

(1.73 ±0.15) X 10i4 

0.084 ± 0.008 

2.85 ±0.20 

1.01 ±0.06 

3.72 ± 0.51 


Note. — Columns list the cluster rsoo, gas mass, total mass derived from the hydrostatic equilibrium equation (Equation (7)), gas 
fraction, spectroscopic-like temperature within (0.15 - 1) X rsoo, and bolometric X-ray luminosities within 245 kpc radii and extrapolated 
to within rsoo. 


TABLE 6 

Physical Properties Derived from the Yy^-M scaling relation 


Subcluster 

rsoo.Yx 

Tfg,500,Yx 

Tf500,Yx 

/gas.Yx 

^^spec 

±X,245 

-('X.rsoo 


(kpc) 

(Me) 

(Me) 


(keV) 

(1043 erg s“l) 

(1043 erg 

G345N 

771 ± 47 

(1.62 ± 0.38) X 10i3 

(1.37 ±0.25) X 10i4 

0.117 ±0.015 

2.08 ±0.17 

1.35 ±0.08 

2.72 ± 0.44 

G345S 

793 ± 42 

(1.33 ± 0.21) X 10i3 

(1.48 ± 0.24) X 10i4 

0.090 ± 0.007 

2.90 ±0.21 

1.01 ±0.06 

3.58 ± 0.48 


Note. — Columns list the cluster rsoOi gas mass, total mass derived from the Yx parameter (Equation (9)), gas fraction, spectroscopic- 
like temperature within (0.15 - 1) X rsoo, and bolometric X-ray luminosities within 245 kpc radii and extrapolated to within rsoo- 




Fig. 5.— Left: Entropy profiles for G345N and G345S. They are remarkably different. G345S presents a very high central 
entropy profile. The corona of the central galaxy has been removed from the profile for G345N. Right: Dimensionless entropy pro¬ 
files for G345N and G345S. We see that the entropy index of G345N is very close to the scaling relation A/iLsoo = 1.42(r/r5oo)^'^ 
at rsoo. On the other hand, G345S exceeds this value at rsoo, suggesting that non-gravitational processes are playing a significant 
role in the thermodynamics of the IGM even at large distances from the center of the cluster. 



where mheated is the gas mass that has been heated (for 
isobaric heating, the heat increment equals the enthalpy 
increment and Ti/Tf = giving the equation 

above). 

A less violent and steady source of energy that could 
explain the high central entropy could be energy injection 
into the ICM from an internal source such as an AGN. We 


can estimate a global entropy by weighting the entropy 
by gas density: 


< K >= 


f Kps dV 

Jp.dV 


(13) 


Within rsoo 800 kpc the ratio of < AT > between 
G345S and G345N is 1.8. From Equation (12), we 
compute that the added heat necessary to increase the 
entropy of G345S is ^ 1.0 x 10®^ erg, which corresponds 
to an AGN power of - 3.3 x lO^'^-lO^® erg s”^ if all this 
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energy has been injected into the ICM within 1-0.1 Gyr, 
respectively. The 843 MHz radio image from the SUMSS 
survey (Bock et al. 1999) of the G345 field shows a bright 
radio source associated with one of the elliptical galaxies 
(ESO 187-lG 025 NED05 in G345S). Its radio power is 
PiA GHz = 1-4 X 10^^ W Hz“^, which corresponds to a 
cavity power of roughly lO'^^-xlO"^® erg s“^ (Birzan et al. 
2008; O’Sullivan et al. 2011). Such an AGN would need 
to inject energy into the ICM for a couple of Gyrs to 
enhance the central gas entropy by a factor of 1 . 8 , which 
also makes this scenario possible. Thus, either an AGN 
sustained over a few Gyrs or the merger of groups could 
increase the central entropy of G345S. 

Departure from the scaling relation K/K^qq = 
1.42(r/r5oo)^'^ (Pratt et al. 2010) is indicative of non- 
gravitational processes, where ATsoo is computed by: 


^500 = 106 keVcm^ (—- 

V 1014 Vo'Mo 



where /b = 0.15 is the baryon fraction (this is the 
assumed value for the baryon fraction in the work of 
Pratt et al. 2010). The right panel of Figure 5 shows the 
dimensionless entropy profile of both G345 subclusters. 
We see that the entropy index of G345N is very close 
to the scaling relation K/K^qq = 1.42(r/r5oo)i i around 
rsoQ. On the other hand, G345S largely exceeds this value 
at rooo, suggesting that non-gravitational processes are 
playing a significant role in the thermodynamics of the 
ICM even at c^qq, supporting the scenario that multi¬ 
ple mergers of groups have boosted the entropy index of 
G345S through shocks. 


6 . PLANCK DETERMINED AND EXPECTED 
SUNYAEV-ZEL’DOVICH SIGNALS 

In this section, we compare the measured and expected 
Sunyaev-Zel’dovich signals from G345. 

The total Sunyaev-Zel’dovich signal is given by the in¬ 
tegral of the Compton parameter, Y = J y dft, where ft 
is the solid angle, and the Compton parameter y is given 
by: 

y = [kTe{r)ne{r) dl, (15) 

where ctt is the Thomson cross section, nieC^ is the elec¬ 
tron rest mass energy, / is the distance along the line 
of sight, and k is the Boltzmann constant. The total 
Sunyaev-Zel’dovich signal also can be expressed as: 


Y = 




(16) 


where Da is the angular size distance of the cluster, 
and P is the electron pressure, P = UekTe. The 
spherical^ Sunyaev-Zel’dovich signal measured by the 
Planck mission is K = 0.0109 ± 0.0032 arcmin^ within 
an angular size of 0 = 5 x d^oo = 118.59 arcmin, 
where 05 oo is the angular size corresponding to rsoo 


^ The Planck Collaboration performs the integral within a 
sphere, instead of performing it along the line of sight to infinity 
(cylindrical integral). 


(Planck Collaboration et al. 2011). 6 corresponds to 6.3 
Mpc at the cluster redshift, which leads to rsoo = 1260 
kpc, and M 500 = 5.94 x IO^^Mq for the entire cluster. 
Based on the X-ray data, we compute rsoo 700-800 
kpc and M 500 1-2 x IO^^Mq for each subcluster. 

Figure 6 shows the P/ancA:-reconstructed Ysz map 
of G345, overlaid with the X-ray isointensity contours. 
Although we see a clear off-set between the X-ray 
and the Sunyaev-Zel’dovich signal peaks, this is con¬ 
sistent with Planck’s much courser spatial resolution 
and higher instrumental noise compared to Chandra. 
Planck Collaboration et al. (2013) showed that an off-set 
as large as 5' between the X-ray and Sunyaev-Zel’dovich 
signal peaks can be expected in the reconstructed j/-map 
for low signihcance objects, due to astrophysical contri¬ 
butions and noise fluctuations. They also showed that 
the SZ signal can be better reconstructed assuming pri¬ 
ors from other wavelengths, such as position, relative in¬ 
tensity between the subclusters, and size. 



Fig. 6 .— P/ancfc-reconstructed Ysz map of G345, over¬ 
laid with the XMM-Newton X-ray isointensity contours. Al¬ 
though we see a clear offset of 10' between the X-ray and 
the Sunyaev-Zel’dovich signal peaks, this is consistent with 
Planck’s much courser spatial resolution and noise, compared 
to Chandra. 

To estimate the expected integrated Compton parame¬ 
ter, we used the electron density and temperature profiles 
we determined for both subclusters to estimate the elec¬ 
tron pressure which we integrate along the line of sight, 
P = where the indices N and S 

correspond to G345N and G345S. The centers of the sub¬ 
clusters are separated by 400 kpc on the sky to match the 
X-ray observations. We computed the integrated Comp¬ 
ton parameter over a sphere of 1260 kpc (the rsoo given 
by the Planck Collaboration) centered between the two 
subclusters (the resulting 2D map is presented in Figure 
7). The Planck Collaboration et al. (2011) determines 
the integrated Compton parameter within Srsoo to be 
Y 5 R 500 — l-SlEsoo, where Ysoo is the integrated Compton 
parameter within rsoo- Using the X-ray derived parame- 
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ters for the gas temperature and density, we computed a 

Y Compton parameter oi Y = 0.0136 ±0.0011 arcmin^, 
which is consistent within ^ O.Scr of the Planck Collab¬ 
oration measured value oiY = 0.0109 ± 0.0032 arcmin^. 

We also modified the pressure model from the univer¬ 
sal profile (Arnaud et al. 2010) that was used to com¬ 
pute Y by the Planck Collaboration to a two compo¬ 
nent model that keeps the relative amplitude of the 
two components fixed to our X-ray model. Using the 
Matched Multi Filters extraction algorithm, we obtained 

Y = 0.0119 ± 0.0025 arcmin^, which is consistent within 
^ 0.6(7 of our X-ray measured value. This value is also 
consistent with the previous Sunyaev-Zel’dovich signal 
measured assuming a universal pressure profile. 



5.87e-09 

5.33e-09 

4.80e-09 

4.27e-09 

3.73e-09 

3.20e-09 

2.67e-09 

2.13e-09 

1.60e-09 

1.07e-09 

5.33e-10 

O.OOe+00 


Fig. 7. — 2.52 Mpc x 2.52 Mpc image of the spherically 

integrated Compton parameter Y with 1 kpc x 1 kpc resolu¬ 
tion. The values are in units of arcmin^. The centers of the 
subclusters were separated by 400 kpc on the sky to match 
the X-ray observations. We see that most of the Sunyaev- 
Zel’dovich signal comes from the northern component, since 
its gas density (and therefore pressure) is considerably higher 
than that of G345S. 


7. A DYNAMICAL MODEL FOR THE G345 
SYSTEM 

In this section we apply the dynamical model presented 
by Beers et al. (1982) to the G345N and G345S system 
to evaluate the dynamical state of the subclusters. 

For the case where the subclusters are gravitationally 
bound, we write the equations of motion in the following 
parametric form; 


R = 

^(1- 
2 ^ 

cosy), 

(17) 

t = 1 

f Ri. \ 

1/2 

(X - siny). 

(18) 

[SGM J 

V = 

f2GM' 

\ smy 

(19) 

V Rm . 

) (1 -cosy)’ 


where Rm is the separation of the subclusters at maxi¬ 
mum expansion, M is the total mass of the system, and 


X is the development angle used to parametrize the equa¬ 
tions. For the case where the subclusters are not gravi¬ 
tationally bound, the parametric equations are: 


^=1/2 (coshy 1), 

^OO 

(20) 

GM, . ^ 

T/3 (sinhx X), 

(21) 

y = w S“hy 

(coshy - 1)’ 

(22) 


where 14o is the asymptotic expansion velocity. The ra¬ 
dial velocity difference, 14 -, and the projected distance, 
Rp, are related to the system parameters by 

Vr = V sin a, Rp = R cos a. (23) 

The total mass of the system is M = (4.03 ± 0.33) x 
10 ^"^ Mq (sum of the masses of both subclusters within 
?’ 20 o)- We assume that the subclusters’ velocities are the 
line of sight velocities of their central dominant galax¬ 
ies. We take Rp = 0.4 Mpc, the projected distance on 
the plane of the sky between the dominant galaxies of 
each subcluster (see Figure 1). The redshift difference 
between these galaxies yields a radial velocity difference 
of Vr = 1134 ± 66 km s“^ (the giant elliptical in G345N 
is ESO 187-G026, at z = 0.047176 ± 0.000143, and in 
G345S the giant elliptical is ESO 187-IG025 NED04, at 
z = 0.043223±0.000180, Smith et al. 2004). We close the 
system of equations by setting t = 12.86 Gyr, the age of 
the Universe at the redshift of these clusters {z = 0.045). 
These equations are then solved via an iterative proce¬ 
dure, which determines the radial velocity difference Id¬ 
as a function of the projection angle a. 

Simple energy considerations specify the limits of the 
bound solutions: 

Vr^Rp < 2GM sin^ a cos a. (24) 

Figure 8 shows the projection angle a as a function of 
the radial velocity difference U between the subclusters. 
The uncertainties in the measured line-of-sight velocity 
and total mass of the system lead to a range in the solu¬ 
tions for the inclination angles (ainf and asup)- The rel¬ 
ative probabilities of these solutions are then computed 
by: 

Pi = COSO; da, (25) 

where the index i represents each solution. Finally, the 
probabilities can be normalized by Pi = Pi/(Yhj^Pi)- 

The radial velocity difference of the subclusters Vr = 
1134 ±66 km s“^ yields two bound solutions and one un¬ 
bound solution for a. For the bound solutions, the sub¬ 
clusters are either approaching each other at 1173 km 
s“^ (17% probability) or at 2636 km s“^ (83% probabil¬ 
ity). The former solution corresponds to an encounter in 
less than 1.2 Gyr, given their separation of ^ 1.56 Mpc. 
The latter corresponds to an encounter in less than 170 
Myr, given their separation of ^ 440 kpc. The unbound 
solution (0.05% probability) corresponds to a separation 
of ^ 17 Mpc. These solutions are presented in Tables 7 
and 8 . Given its low probability, the unbound solution 
can be neglected, while the bound solution in which the 
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Fig. 8 .— Projection angle a as a function of the radial velocity 
difference Vr between the subclusters. BO, BI, and UO stand for 
Bound Outgoing, Bound Incoming, and Unbound Outgoing solu¬ 
tions, respectively. Solid blue and red lines correspond to bound 
and unbound solutions, respectively. The vertical solid line corre¬ 
sponds to the radial velocity difference between the giant ellipticals. 
Dashed lines correspond to 68% confidence ranges. 

separation between the clusters is ^ 440 kpc is highly 
favored (83% probability). 

We note, that the dynamical analysis method from 
Beers et al. (1982) assumes a purely radial infall (no 
angular momentum). Also, the way the probabilities 
are computed favors small angle solutions. The small 
angle bound solution (25.48°) is highly favored {P = 
83%) compared to the other bound solution (75.13° 
- P = 17%), despite its supersonic infalling velocity 
(^ 2600 km s“^ - the sound speed of a 3 keV cluster 
is ^ 900 km s“^, therefore 2600 km s“^ corresponds to a 
Mach number of ^ 3). The virial radius of these subclus¬ 
ters is roughly 1.2 Mpc. If, indeed, they were separated 
by only ^ 440 kpc, moving at ^ 2600 km s“^, shock 
discontinuities in the X-ray surface brightness would be 
seen in the region between their cores. Furthermore, fit¬ 
ting the spectrum in a rectangular region (105 kpc x 
345 kpc) between the subclusters gives kT = 3.30to;53 
at 68% confidence, providing no evidence of shock heated 
gas between the subclusters. In this analysis, the prob¬ 
abilities for the bound solutions should be treated with 
caution, as we have no information about the angular 
momentum of this system. 


TABLE 7 

Best Fit Parameters for the Bound Incoming Solutions of 
THE Dynamical Model. 


X 

a 

R 

Rm 

F 

P 

(rad) 

(degrees) 

(kpc) 

(kpc) 

(km s“^) 

(%) 

4.952 

75.13 

1559.5 

4089.6 

1173.1 

17 

5.600 

25.48 

443.1 

3952.9 

2635.6 

83 


Note. — Columns list best fit for the x S'Ud a for the bound 
solutions of the dynamical model, and the corresponding values for 
it, itm, F, and probability of each solution. 


TABLE 8 

Best Fit Parameters for the Unbound Outgoing Solution 
OF THE Dynamical Model. 


X 

(rad) 

a 

(degrees) 

R 

(kpc) 

U Uoo 

(km s“^) (km s“^) 

P 

(%) 

3.157 

88.67 

17213.7 

1134.1 1041.5 

0.05 

Note. 

— Columns list best fit for the x ^ for the unbound 


solution of the dynamical model, and the corresponding values for 
it, V, Foo, and probability of this solution. 

7.1. A Modified Dynamical Model 

Since the Beers et al. (1982) dynamical model is based 
on the timing argument (see e.g. Kahn & Woltjer 1959), 
which Li & White (2008) have shown to be biased and 
over-constrained (failing to reproduce the scatter ob¬ 
served in N-body simulations), we also investigate G345 
using a modified version of the Dawson (2013) Monte 
Garlo dynamic analysis method®. The Dawson (2013) 
method relaxes many of the constraints in the timing 
argument and examines all merger scenarios consistent 
with the observed state, enabling it to capture the ob¬ 
served scatter in the N-body simulations. Additionally 
this method models the two subclusters as NFW halos 
(Navarro et al. 1996) so we can more accurately® esti¬ 
mate the time-till-collision (TTC), the period between 
collisions (T), and the eventual relative collision veloc¬ 
ity (V 3 D(tcoi)) (for a precise dehnition of these quan¬ 
tities, please refer to Dawson 2013). We modified the 
Dawson (2013) method slightly, since it is designed for 
post-merger systems and there is strong evidence that 
G345 is a pre-merger system. We removed the prior that 
the time-since-collision^ (TSC) be less than the age of 
the Universe, since TSC is now recast as the TTC. We 
also now allow for unbound scenarios, see Equation (24), 
however we require that 14 be less than the line-of-sight 
Hubble flow velocity for the unbound realization to be 
considered valid, 

Hubble — H{z)Rsin{a), (26) 

where H (z) is the Hubble parameter at the average red- 
shift of the two subclusters. We summarize the results 
of this analysis in Table 9. 

From Figure 9 it can be seen that the large rela¬ 
tive velocity and small separation highly favored by the 

® Monte Carlo Merger Analysis Code (MCMAC Dawson 2014). 

® To within ~5-10% agreement with N-body simulations 
(Dawson 2013). 

^ We define the time of collision to be the time of the first peri¬ 
centric passage. 
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Beers et al. (1982) dynamic model are disfavored by the 
modified Dawson (2013) dynamic model, as well as being 
inconsistent with the observed gas properties. However, 
there are still a number of realizations in the Monte Carlo 
analysis where effects of the merger on the gas might 
be expected. In principle these also could be excluded 
from the posterior distributions. Based on this analy¬ 
sis, we find that the subclusters are likely to collide in 
0.5 ± 0.2 Gyr with a relative collision velocity of 2000 ± 
100kms“^, see Figure 10. Based on this dynamic analy¬ 
sis, there is a slightly larger probability (2.6% vs. 0.05%) 
that the subclusters are unbound, but this scenario is 
still unlikely. The unbound parameter estimates are also 
summarized in Table 9. 



y{tohs) (km s ^) 


Fig. 9.— The marginal posterior of the three-dimensional sub¬ 
cluster separation and relative velocity ^30 in the observed 
state tobs? inferred from the modified Dawson (2013) dynamics 
analysis. Dark and light blue contours represent the 68 % and 95% 
confidence regions, respectively. 

8 . CONCLUSIONS 

We have presented Chandra, XMM-Newton, and 
ROSAT observations of G345, a Planck detected dou¬ 
ble cluster, and provided measurements of temperature, 
gas and total masses, gas fraction, entropy profiles, ex¬ 
pected Sunyaev-Zel’dovich signal, and X-ray luminosi¬ 
ties. Both the north and south subclusters have X-ray 
surface brightnesses that are well described by a /3-model 
within rsoo, with /3 ~ 0.5. Both subclusters have gas 
masses in the range 1-2 xIO^^Mq and total masses in 
the range 1-2 x IO^'^Mq, and gas mass fractions in agree¬ 
ment, within the confidence range, of those found by 
Vikhlinin et al. (2006) for clusters with similar total mass 
(0.12±0.02 for the northern subcluster and 0.09±0.01 for 
the southern one). We show that the G345 subclusters 
that are very likely (> 97% probability) gravitationally 
bound and infalling and will collide in 500 ± 200 Gyr. 
We show that there is good (0.8cr) agreement between 
the expected Sunyaev-Zel’dovich signal predicted from 
the X-ray data and the measured value from the Planck 


1.4 


1.2 



0.2 


1600 1800 2000 2200 2400 

KAoi) (km s“^) 

Fig. 10.— The marginal posterior of the relative three- 
dimensional subcluster velocity at the time of collision 
and the time it will take to collide given their observed state, in¬ 
ferred from the modified Dawson (2013) dynamics analysis. Dark 
and light blue contours represent the 68 % and 95% confidence re¬ 
gions, respectively. 

mission, and 0.6cr agreement when the Planck value is 
re-computed assuming a two component pressure model, 
with relative amplitudes fixed using the X-ray results. 
The high central entropy in G345S can be explained ei¬ 
ther by the mergers of several groups, as suggested by the 
presence of five massive elliptical galaxies or by AGN en¬ 
ergy injection, as suggested by the presence of a bright 
radio source in the massive elliptical galaxy ESO 187-IG 
025 NED05, or by a combination of both processes. 
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TABLE 9 

Parameter Estimates from Modified Dawson (2013) Dynamics Analysis 


Parameter 

Units 

1 

Bound Scenarios {C= 

=97.4%) 

1 

Unbound Scenarios {C= 

=2.6%) 




Location^ 

68 % LCL-UCLb 

95% LCL-UCLb 

Location®’ 

68 % LCL-UCLb 

95% LCL-UCLb 

M200,N 

1014 Mq 

1.6 

1.4 - 1.8 

1.1 - 2.0 

1.6 

1.3 - 1.8 

1.1 

- 2.0 

^^200,3 

IQi^M© 

2.5 

2.3 - 2.8 

2.0 - 3.0 

2.5 

2.3 - 2.7 

2.0 

- 3.0 

2 N 


0.0472 

0.0470 - 0.0473 

0.0469 - 0.0474 

0.0472 

0.0470 - 0.0473 

0.0469 

- 0.0475 



0.0433 

0.0431 - 0.0434 

0.0429 - 0.0436 

0.0432 

0.0430 - 0.0434 

0.0429 

- 0.0436 

Rp 

Mpc 

0.40 

0.36 - 0.44 

0.32 - 0.48 

0.41 

0.37 - 0.45 

0.33 

- 0.50 

a 

degree 

59 

46 - 70 

37 - 77 

89 

89 - 90 

89 

- 90 

R 

Mpc 

0.81 

0.60 -1.3 

0.48 -2.0 

30 

20 - 70 

15 

- 160 

Rm 

Mpc 

2.0 

1.6 - 2.9 

1.4 -5.7 





^D(^obs) 

km 

1300 

1200 - 1600 

1100 - 1900 

1100 

1100 - 1200 

1000 

- 1300 

Fid (Col) 

km 

2000 

1900 - 2100 

1800 - 2300 





TTC 

Gyr 

0.5 

0.3 - 0.7 

0.2 - 1.1 





T 

Gyr 

4.9 

3.7 -7.8 

3.1 - 15 






® Biweight-statistic location (see e.g. Beers et al. 1990). 

Bias-corrected lower and upper confidence limits, LCL and UCL respectively (see e.g. Beers et al. 1990). 
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APPENDIX 

STELLAR AND LOW MASS X-RAY BINARY CONTRIBUTIONS TO THE X-RAY LUMINOSITY OF 

ELLIPTICAL GALAXIES 

Stellar Contribution 

Revnivtsev et al. (2007) showed that the unresolved X-ray halo in M32 can be best explained by the cumulative 
emission from cataclysmic variables and coronally active stars. In a following work, Revnivtsev et al. (2008) used deep 
Chandra observations to measure the unresolved X-ray emission in the elliptical galaxy NGC 3379. They suggested 
that the old stellar populations in all galaxies can be described by a universal value of X-ray emissivity per unit stellar 
mass or K-band luminosity. 

From the 2MASS K-band image of the giant elliptical in G345N (ESO 187-G026), we extract its K-band luminosity 
and determine its stellar mass. The K-band luminosity measured within 7 kpc from the center of ESO 187-G026 (for 
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comparison with the galaxy’s X-ray luminosity measured within the same region) is Lk = 1-6 x 10^^ erg s Bell et al. 
(2003) showed that the relation between the K-band luminosity and the stellar mass is given by: 

log = OK + 6 k X (color) (Al) 

where the M/L ratio is given in solar units. For the B — R color, ok = -0.264 and 6 k = 0.138 (Bell et al. 2003). We 
calculate B — R = 1.63 for ESO 187-G026, which gives M = 1.78 x IO^^Mq. 

The X-ray luminosity in the soft (0.5 - 2.0 keV) band due to stellar emission from ESO 187-G026 can now be 
computed using the relation given by Revnivtsev et al. (2007): 


Aq. 5-2.0 keV — 7 X 10' 


,38 


/ M 


V 1 O 11 M 0 


erg s 


(A2) 


Using Equation (A2) we compute a soft X-ray luminosity of Ao. 5 - 2.0 keV = 1-25 x 10^® erg s 0.5 % of the total 
X-ray emission of the galaxy. 


Low Mass X-ray Binaries 

Low mass X-ray binaries (LMXB) also contribute to the unresolved X-ray emission of galaxies. Gilfanov (2004) 
showed that the distribution of near-infrared light in all galaxies closely traces the azimuthally averaged spatial 
distribution of LXMBs. To describe it quantitatively, they defined a template for the X-ray luminosity function as a 
power law with two breaks: 


r Ki (Aas/Lb^i) , Lag < L;,p 

= J AI 2 (A38/L{,_2) , A{,_i < Lgg < Lf,_2 

dAag I LTg (Lgg/L 

cut ) ^ j A{,_2 < Agg < Lent 

0; Agg > Acut 

where Lgg = L^/IO^® erg s“^ and normalizations Ki,K 2 , and ATg are related by 

K 2 = K, (L5,i/L5,2)“^ 

ATg = K 2 (Ab,2/Acut)°'^ 

They fixed the value of the high luminosity cut-off at Lcut = 500 x 10^® erg s“^. Due to the steep slope of the luminosity 
function above A{,_ 2 , the results are insensitive to the actual value of Lcut- Typically, for nearby galaxies, the source 
detection threshold defines the luminosity cut-off. However, for ESO 187-G026 we cannot resolve any point sources, 
so the luminosity cut-off is set to their fixed maximum value of Lcut = 500 x 10®® erg s“®. 

Gilfanov (2004) provides the best fit value for the average normalization of ATg = 440.4 per 10®® Mq, and the best 

fits to the parameters of the luminosity distribution of ag = 1.0, 0:2 = 1.86, ag = 4.8, Lb,g = 0.19, and Lb ,2 = 5.0. 

We can now compute the cumulative number of sources using: 

iVx(> L) = / ^ dLgg, (A5) 

J L U-Ogg 

and the total luminosity by: 

dN 

Ax(> A) = / —— Lgg dLgg, (A6) 

JL Utvgg 

which yields Lq. 5 - 2.0 kev(> 10^^ erg s“®) = 1.62 x 10®° erg s“® for ESO 187-G026. This leads to a (stellar -|- LMXB) 
luminosity of Lq, 5 - 2.0 keV = 1-74 x 10®° erg s“®, which corresponds to ~ 7% of the galaxy’s total 0.5-2.0 keV measured 
X-ray luminosity within 7 kpc of its center. 


(A3) 


(A4) 






